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The Foldy-Wouthuysen transformation is known to uncover the nonrelativistic limit of a gen¬ 
eralized Dirac Hamiltonian, lending an intuitive physical interpretation to the effective operators 
within Schrodinger-Pauli theory. We here discuss the opposite, ultrarelativistic limit which requires 
the use of a fundamentally different expansion where the leading kinetic term in the Dirac equa¬ 
tion is perturbed by the mass of the particle and other interaction (potential) terms, rather than 
vice versa. The ultrarelativistic decoupling transformation is applied to free Dirac particles (in the 
Weyl basis) and to high-energy tachyons, which are faster-than-light particles described by a fully 
Lorentz-covariant equation. The effective gravitational interactions are found. For tachyons, the 
dominant gravitational interaction term in the high-energy limit is shown to be attractive, and equal 
to the leading term for subluminal Dirac particles (tardyons) in the high-energy limit. 
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I. INTRODUCTION 

The Foldy-Wouthuysen transformation [l| is an es¬ 
tablished method, used to calculate the nonrelativistic 
limit of effective Hamiltonians describing spin-1/2 par¬ 
ticles. The procedure has been applied with good effect 
to the Dirac-Coulomb Hamiltonian ii , uncovering the 
fine-structure terms for atomic levels, notably, the zitter- 
bewegung term, and the Russell-Saunders (spin-orbit) 
coupling (Thomas precession). Recently, the analogue 
of the Russell-Saunders coupling in a gravitational field 
(the Fokker precession, see Refs. (3^) has been recov¬ 
ered from the gravitationally coupled Dirac Hamiltonian, 
which is referred to as the Dirac-Schwarzschild Hamilto¬ 
nian Q. 

The Foldy-Wouthuysen program, in its original 
form [Ij, is inherently perturbative in nature. In a typical 
case, the structure of a generalized Dirac Hamiltonian is 
given as a. ■ p + [3 m + SH (in the standard Dirac repre¬ 
sentation of the Clifford algebra, see Appendix El). Here, 
the “dominant” term is taken as (3 m, where (3 is the 
4x4 Dirac (3 matrix, a ■ p \s the kinetic operator (d? is 
the vector of Dirac a matrices, and p is the momentum 
operator), and 5H contains the potential terms. One 
then expands about a Dirac particle “at rest”, with the 
dominant term given by the “rest mass” term f3m. The 
Foldy-Wouthuysen procedure then uncovers the leading 
nonrelativistic kinetic term a ■ p ^ p'^/(2m) -|- ... and 
transforms the potential terms SH into a form where 
the operators acquire an intuitive physical interpretation. 
At some risk to oversimplification, one can say that the 
Foldy-Wouthuysen transformation applies to the regime 
\a-^^ \l3m\, and \6H\ <C |/3m|. 

In some cases, such as for a free Dirac particle [3, 
it is sometimes possible to perform the transformation 
without any perturbative expansion in the momenta or 
other expansion parameters. There have been attempts 
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to generalize the idea of a nonperturbative method to 
more general Hamiltonians, and a set of interesting iden¬ 
tities have been derived Ref. Q. However, the alterna¬ 
tive Foldy-Wouthuysen transformation Q suffers from 
an explicit breaking of the parity symmetry in the trans¬ 
formation, which involves the fifth current, and is known 
to produce spurious parity-breaking terms in a number 
of applications, e.g., to the Dirac-Coulomb Hamiltonian 
(for an overview, see Refs. Sii). In general, nonper¬ 
turbative methods (in the momenta of the Dirac parti¬ 
cles) can only be applied when considerable additional 
information is available for a specific Hamiltonian under 
investigation, and when additional approximations are 
made, such as the neglecting terms of second order and 
higher in the field strengths [see Eq. (21) of Ref. mi 

To the best of our knowledge, the opposite perturba¬ 
tive expansion, namely, perturbation theory of a Dirac 
Hamiltonian about the ultrarelativistic limit, has not yet 
been considered in the literature; it is the subject of the 
current paper. This expansion has to follow a funda¬ 
mentally different paradigm; in the ultrarelativistic limit, 
mass terms and potential terms are suppressed in com¬ 
parison to the kinetic term; the expansion is valid in the 
regime |/3m| <C \a ■ p\, and \6H\ -C \d ■ p\. Ultrarela¬ 
tivistic particles are best described in the helicity basis 
[Chap. 23 of Ref. [HI] , while in fact, the solutions of the 
free Dirac equation approach those of the Weyl equa¬ 
tion in the massless limit (see Chap. 2.4.3 on p. 87 of 
Ref. i). The Weyl equation describes massless spin-1/2 
particles, which transform under the fundamental (^,0) 
representation of the Lorentz group and travel exactly at 
the speed of light (these are the “neutrinos in the original 
standard model”). 

We here investigate the ultrarelativistic decoupling 
transformation with a special emphasis on the gravita¬ 
tional coupling of a particle to a central gravitational 
field. To this end, in Sec. m we briefly recall the un¬ 
derlying covariant formalism, distinguishing the case of 
a “normal” (subluminal) Dirac particle from a particle 
described by the tachyonic Dirac equation [3|6|- The 
latter equation describes faster-than-light particles, still 
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in a fully Lorentz-covariant formalism [l^ . The ultra- 
relativistic limit specifically is relevant to tachyons be¬ 
cause these particles cannot travel slower than light; they 
remain superluminal upon Lorentz transformation |13l - 
[i3- In the ultrarelativistic limit, the particle’s speed ap¬ 
proaches the light cone and the influence of tardyonic as 
well as tachyonic mass terms are suppressed in compari¬ 
son to the kinetic terms. In Sec. imi the ultrarelativistic 
decoupling transformation is applied to gravitationally 
coupled tardyonic and tachyonic particles. Conclusions 
are reserved for Sec. Hvl 


an ultrarelativistic decoupling (ultrarelativistic Foldy- 
Wouthuysen) transformation. To this end we define the 
energy operator 

£ = -^-p, (6) 

and the transformation (see Sec 4.2 of Ref. 0) and 
Sec. 2.2.4 of Ref. 0) 

5fd = -i /3 7 " ^0 , = 5fd . (7) 

so that the unitary transformation C/pD becomes 


II. FREE PARTICLES 

A. Free Tardyonic Transformation 

In principle it is well known that the Weyl equa¬ 
tion, which describes a massless spin-1/2 particle, splits 
into two equations, describing a left-handed and a right- 
handed spinor (see Chap. 23 of Ref. 0 and p. 87 of 

Ref. 0 ), 


idttpL = Hl'^l , HL = -a-p, ( 1 ) 

idtipR = Hri/jr, HR = a-p, (2) 

The Weyl equations break parity; a left-handed spinor 
transforms into a right-handed solution under the parity 
operation. However, it is well known that the Dirac equa¬ 
tion, whose bispinor solutions are constructed by stack¬ 
ing the helicity spinors on top of each other, conserves 
parity [H]. 

The massless equation, in turn, corresponds to the ul¬ 
trarelativistic limit for a massive Dirac particle; we would 
thus expect that the Dirac equation splits into two equa¬ 
tions, describing left- and right-handed Weyl spinors, in 
this limit. Thus, if we are to recover the massless (Weyl) 
limit, plus corrections, from the Dirac equation, then we 
need to necessarily invoke a parity-breaking transforma¬ 
tion. We start from the free Dirac Hamiltonian 

Hfd = a ■ p + (3m = f ^ ^ (3) 

\(T ■ p —m J ^ ^ 


CpD = exp(i5'FD) = cos 


Choosing 0 so that 


cos I 2 


^0 


sin I 2 


m 0 



.,5 1^ 


/3 7 ^ sin 



1^ 

m 

a/p 2 -I- ’ 


(9a) 

(9b) 


one finally obtains 

"Hfd = C7fd Hfd = — \/p^ +m’^ . (10) 

In explicit (2 x 2)-matrix subform, 

Hfd = - (3 \/p 2 _|_ ^2 

f -Vp^ + m'^ 0 

= 1 ^ ^ ^ 

V ^ + m 2 

it becomes clear that the separation into a left-handed 
(upper diagonal) and a right-handed (lower diagonal) 
Hamiltonian has been achieved. 

The eigenstates of the Hamiltonian (fTOll fulfill 



and invoke the following unitary, parity-breaking trans¬ 
formation 


( 4 ) 

which transforms Hfd into Hfd = U Hfd U ^, 

HFD = -/3S-P + 7^m=(-^-^%"^.) (5) 


The initial rotation with the U matrix corresponds to 
a change of the basis of the Dirac algebra, into the so- 
called Weyl basis (see Appendix!^. The mass terms are 
now off-diagonal and we may try to eliminate them by 


i C) = Hfd V’/C 7/! i = 1,2,3,4. (I2a) 

The first two solutions can be written as 

V'i(t,r)= Q-'Et+itr ^ Hfd tpi = E tpi, 

(12b) 

tp2{t,r)= Q^Et-ik-r ^ HfD'02 = -HV'2 • 

(12c) 

The physical momentum is k, and the helicity eigenvalue 
is negative for both solutions, S • k'ipi ^2 = ~l/’i. 2 , and 
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E = VP + m? (here, k is the unit vector in the k di¬ 
rection). The solution V 2 describes an antiparticle. The 
two solutions of right-handed helicity are 

Mt, r) = ) e-'^ , iJpD = 

(12d) 

Mt, f)= ) e'^ , ifpD V 4 = -£; V 4 , 

(12e) 


The helicity is positive for these two solutions, E- fc ^ 3,4 = 
'03,4; with ipi describing an antiparticle. The eigenvalues 
of the £ operator for Vi, 2 , 3,4 are E,—E,—E,E, respec¬ 
tively. If we apply the formalism to a Dirac neutrino, 
them ipi would describe a left-handed neutrino, V 2 would 
describes a left-handed antineutrino, whereas Vs and V 4 
would describe right-handed neutrinos and right-handed 
antineutrinos, respectively. For completeness, we recall 
the form of the helicity spinors Q, 


a+{k) = 
a-{k) = 


pi V 


cos (I) 

pin (I) • 

' - sin (I) e’ 
cos (I) 


-iip 


(13a) 

(13b) 


where 6 and ip are the polar and azimuthal angles of the 
k vector. 


B. Free Tachyonic Transformation 


astrophysical observations have recently been recorded 
in Refs. [25l - [^ : other theoretical studies concern Dirac 
equations with Lorentz-violating terms which can lead 
to super luminal propagation for neutrinos [ 2 ^, . The 

tachyonic Dirac Hamiltonian has recently been identi¬ 
fied as a pseudo-Hermitian (“y^-Hermitian”) Hamilto¬ 
nian in Ref. [ 2 ^. Independent of the phenomenological 
relevance of the concept of tachyons, the current section 
of our paper, and Sec. IHI HI demonstrate that it is pos¬ 
sible to uncover the leading terms of generalized pseudo- 
Hermitian [sTI - ISSIl Dirac Hamiltonians in the ultrarela- 
tivistic limit using the relativistic decoupling transfor¬ 
mation. 

The accepted generalized Dirac Hamiltonian for a free 
tachyonic Dirac particle is given as [illiiliS 

FItd = a • P +/37^ m, (14) 

which is 7 ®-Hermitian, i/xD = 7 ^^^td 7^- then fol¬ 
low the same procedure outlined in Sec. Ill Al and begin 
by performing the initial rotation U [see Eq. 0], giving 
us 


Htd = VFItd V-i =/3£:-f/ 37 '^m. (15) 

The Hamiltonian Htd is /3-Hermitian, i.e., Htd = 
j3. Here, (3 is the Dirac /3 matrix which takes the 
role of the 7 ® matrix in the Weyl representation (see Ap¬ 
pendix]^. The /3-Hermitian operator S'td in this case 
reads as 

^TD = -i/3^^0 = -i7'^f0, 5 td = /3 5+d/3. 

(16) 

The transformation 


As we have just shown, one may accomplish an ex¬ 
act diagonalization (in spinor space) of the free Dirac 
Hamiltonian using the ultrarelativistic decoupling trans¬ 
formation. However, one might counter argue that this 
result is in principle familiar: An exact diagonalization 
can also be accomplished using the Foldy-Wouthuysen 
transformation (in its original form) for the free Dirac 
Hamiltonian (see Sec. 4.2 of Ref i). The ultrarelativistic 
transform leads to a form which asymptotically is equal 
to the Weyl Hamiltonian (helicity basis) of a massless 
particle, as it should be (in the ultrarelativistic limit). 
Here, we shall make the point that, unlike the original 
Foldy-Wouthuysen transform, which can only be applied 
to tardyons, the ultrarelativistic decoupling can also be 
used for Lorentz-invariant tachyons [12, Oj 113 j whose 
velocity remains superluminal upon Lorentz transforma¬ 
tion [l3-[i3. 

A few general remarks on tachyonic spin-1/2 particles 
might be in order. The tachyonic neutrino hypothesis 
remains one of the driving forces behind the study of 
the tachyonic Dirac equation [I2|. The algbraic struc¬ 
tures underlying the tachyonic spin-1/2 equation have 
recently been studied in greater depth (see Refs. [20l - [^ 
and references therein). Pertinent potentially relevant 


Ftd = exp(iS'TD) = cosh 





( 17 ) 


fulfills the identity 

P Utd = exp(iS'TD)/3 exp(iS'TD) = P , (18) 


i.e., it is /3-unitary. It therefore conserves the VT- 
symmetric scalar product (V|/3|V)- Choosing [redefining, 
see Eq. 0] 0 so that 


one obtains 
■Htd = 


cosh j 

{^me\ 


(19a) 

K b1 J 

— m? 

sinh j 


m 

(19b) 

\ b1 ) 

y/jP — np 


C/td Htd = /3—. (20) 

IpI 


This amounts to the exact ultrarelativistic decoupling 
transformation of the free tachyonic Hamiltonian, in the 
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helicity (“Weyl”) basis which has been shown to lead 
to a very efficient description of the tachyonic bispinor 
solutions [H [ 33 . The Taylor series expansion of 
"Htd gives rise to the terms 


Htd 




= -/3S.p 





( 21 ) 


which are the correction terms in the high-energy limit. 
For the tachyonic case, the eigenstates of the Hamilto¬ 
nian (1^01) are still given by Eq. (HH), but one has to re¬ 
place E = \/k‘^ + m? —>■ in the sense of the 

tachyonic dispersion relation. 


III. TRANSFORMATION WITH 
GRAVITATIONAL COUPLING 


A. Gravitational Tardyonic Transformation 


Schwarzschild radius Vg (first order in G) and keep in¬ 
verse powers of £ up to order \/£. 

It is advantageous to write the Hamiltonian (l2^ as 

HDS = /3f-f {f,^}+ODS. (26) 

The first commutator is given as 


[^ds,Hds]= [^DS,/3f]- 

+ [S'ds, C’ds] • 


{^^ 7 } 


(27) 


Let us investigate the first commutator for 

which one finds after a somewhat tedious calculation, 


[-S'ds,/?^] = IOds + [£, [^,C>Ds]] ^ , (28) 

where the double commutator is proportional to a three- 
dimensional Dirac-(5 function plus a spin orbit coupling 
term. 


The study of the gravitationally coupled Dirac equa¬ 
tion, for massless particles, was initiated by the ques¬ 
tion of how the neutrinos (assumed by symmetry to 
be strictly massless in the Original Standard Model) 
interact with gravitational fields [s^- We follow this 
route and start from the gravitationally coupled Dirac- 
Schwarzschild Hamiltonian [3 

i?DS = ^{a-71-^}+^m(l-^) • (22) 


After the initial transformation into the Weyl basis, one 
finds for Hds = U Hds [see Eq. (|4])] 



In order to proceed with the ultrarelativistic decoupling 
transformation, we identify the odd part Ods of Hds and 
define 


Ods = 7 ® m ^I - , S'ds = -i^ |c>DS, , (24) 

where 5 'ds is Hermitian. The unitary transforma¬ 
tion C/ds = exp(iS'Ds) is applied to calculate T-Lds = 
Utis Hds perturbatively, 

Hds ~ Hds + [“^ds , Hds] + ^ [-S'ds , [^ds , Hds]] + ■ • ■ 

(25) 

which is a series of nested commutators, as with the clas¬ 
sic Eoldy-Wouthuysen transformation [l|. In the fol¬ 
lowing, we carry the calculation to first order in the 


[f, [f, Ods]] = - 271 r, 5 ( 3 ) - r, ^ , (29) 

which is of order unity in the expansion in inverse powers 
of £. Despite the fact that the double commutator has 
two instances of the operator £, the commutators ensure 
that these instances of £ act only on OdSi and not on 
the reference state wave function, which would otherwise 
generate inverse powers of £. Eor the operator £ (or the 
inverse thereof) to be the “dominant term”, it must oper¬ 
ate on a wave function describing a high-energy particle. 
Thus 

i[£,[£,Ous]]^ = o(J^^^0. (30) 

Alternatively, one may observe that, when using the Weyl 
free-spinors given in Eq. (USD as reference states, the the 
expectation values of both the Dirac-5 function and the 
spin-orbit coupling term (S-L/r^) vanish for both diago¬ 
nal as well as off-diagonal matrix elements. In conclusion, 
to the order relevant for our investigation, we can replace 

[S'ds,/?^’] IOds , (31) 

in our approximation. This relation ensures the odd 
terms will be canceled out to the first order in Ods when 
calculating the transformed Hamiltonian 'Hds according 
to Eq. (1^51) . One also establishes that 



(32b) 
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so that the first commutator becomes 


B. Gravitational Tachyonic Transformation 


['S'ds,Hds] = i ( C’ds + — 


The double commutator is then 


m? P rrp fir. 


E' r 
(33) 


[‘S'ds, [5'ds,Hds]] = i ( [<S'ds,C>ds] + 5'ds,7^™—^ 


5'ds,/3 


m 


21 


5, 


P w? 


i II 

£ ’ r 


(34) 


where the first term is known from Eq. (I32bl) . The other 
relevant commutators are 


c 5 ■o'^ / 1 i. 

<5ds, 7 "1— i f > 


£ r 


Sus.p"^ 


= 01^1^0, 


S'ds, 2 


We then have 


i H 

£’ r 


= o 


0 . 


[Sds, [S, Hds]] = p— - Pm^ \ 


(35a) 

(35b) 

(35c) 

(36) 


where again we neglect higher-order terms. Because Sds 
carries an inverse power of we can neglect the triple 
commutator, 


[^Ds, [5 'ds, [*S'ds, Hds]]] = O 


0 . 


(37) 


Thus 


Hds = Hds + i [•S'ds , Hds] + ^ [<S'ds , [S'ds , Hds]] 

where 

Ods = -7'^™-- (39) 

r 

The second iteration of the transform with 

‘S'ds = ; t^DS = exp(i5'Ds)) (40) 

will serve only to eliminate the remaining odd term. 
Thus, the final result for Hds = C^ds ^iis ^ds^ reads 
as 

,4.) 

for a gravitationally coupled high-energy Dirac particle. 


We start from the tachyonic, gravitationally coupled 
(TG) Dirac Hamiltonian derived in Appendix [C] 

H^G = \[a-p,(l-'^)]+Pl^m{l-^Y ( 42 ) 

The initial rotation into the Weyl basis of the Dirac alge¬ 
bra using the transformation U defined in Eq. leads 
to the Hamiltonian Htg = U Htg U~^, which reads as 


HTG = f {£,l-(i} + /37">n(l-^) 
\ -"»(!-if) 


(43) 


where £ has been defined in Eq. One identifies the 
odd part of the Hamiltonian Htg £md writes 


Otg = Pl^ mil -, S'tg =|c>TG, • 

. . 

The /3-unitary transformation C/tg = exp(iS'TG) is ap¬ 
plied to calculate HJpQ = C/tg Htg C^tG’ perturbatively, 


Htg = Htg + ['^'tg, Htg] + ^ [<S'tg, [-S'tg , Htg]] + ■ • ■ 

(45) 

in full analogy with the Dirac-Schwarzschild Hamilto¬ 
nian. After a somewhat tedious calculation, neglecting 
(as before) the Dirac-(5 and spin-orbit coupling terms, 
one finds 


['S'tg, P£] = iC^TG , 




= i P 7^m ■ 


m 


[S'tGjC’tg] = ^ ~ i/3 2 

The first commutator becomes 


i If. 

f ’ r 


(46a) 

(46b) 

(46c) 


[S'tg, Htg] = i C^tg — P'y^m— + P— 

r t 


Pm? 


i !a 

£' r 


(47) 


The double nested commutator is 
[^tg, ['S'tg, H]] = i (['S'tg, C>tg] - 


5tg,/3^ 


STG,PTrn- 


1 r. 


STG,P—i^,f 


, (48) 


The last two commutators are of order l/£^ and can 
therefore be neglected. With the help of the result 


S', PY m — 

r J 


, PrrP 


i a 

£’ r 


( 49 ) 
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and with Eq. (I46cl) . one finds 
[S'tg, [5'tg,Htg]] = +/3w^ 1^, • (50) 

Thus, 


Htg = Htg + i ['S'tG: Htg] + ^ [-S'tG) [5'tGj Htg]] 

where 

O' = m— . (52) 

r 

A second transformation with 

‘5'tg = |c>xgj ! t^TG = exp(iS'xG) > (53) 


eliminates OxG’ obtain the following final result 

for i^TG = U!^g Htg tg> 


Htg 



(54) 


It differs from the result given in Eq. (I4ip only in the 
sign of the kinetic term —m^/( 2 £l), due to the tachyonic 
dispersion relation. 


IV. CONCLUSIONS 

We have studied the ultrarelativistic decoupling trans¬ 
formation for the free Dirac equation (Sec. Ill AT) , and for 
the free tachyonic Dirac equation fSec. lIIBll . These trans¬ 
formations lead to a full separation of the Dirac equation 
in the helicity basis. Unlike the exact Foldy-Wouthuysen 
transformation, which transforms the free Dirac Hamil¬ 
tonian into the form /3 \/j? + rm? (see Ref. i) , the ul¬ 
trarelativistic transformation leads to a separation in the 
helicity basis, with the transformed Hamiltonian being 
proportional to (—/3S [see Eqs. (fTOl) and ([20])]. The 

eigenstates of this Hamiltonian are naturally obtained in 
the helicity basis [see Eq. (in])] and are formally iden¬ 
tical (upon a redefinition of the energy parameter E) 
to the eigenstates of the massless Dirac equation (see 
Chap. 2.4.3 on p. 87 of Ref. Q). The latter eigenstates 
are known to transform under the fundamental (^, 0 ) 
representation of the Lorentz group; the “helicity of the 
massless spinors does not flip upon a Lorentz transforma¬ 
tion”. This observation is intimately linked to the fact 
that massless Dirac spinors describe particles which al¬ 
ways move at the speed of light; it is impossible to “over¬ 
take” the particle, which otherwise leads to a helicity flip 
(see Ref. |^b 

The initial unitary transformation U given in Eq. o 
transforms the Dirac equation into the Weyl basis (see 


Appendix , which is naturally identified as the ultra¬ 
relativistic basis for the description of the Dirac algebra: 
Namely, the Dirac o? matrices are replace, in the Weyl 
basis, by matrices (—/3E • p), which are diagonal in the 
(2 X 2)-spinor space and describe the Hamiltonian for a 
massless Dirac particle. The Dirac and Weyl representa¬ 
tions of the Clifford algebra are complementary: In the 
Dirac basis, the “dominant term” in the Hamiltonian is 
/3m, and the odd (off-diagonal) kinetic terms a-p are 
eliminated by the Foldy-Wouthuysen transformation. In 
the Weyl basis, the kinetic term {—/3Yi-p) is diagonal 
(“dominates in the ultrarelativistic limit), while the off- 
diagonal mass term 7 ® m needs to be eliminated by the 
ultrarelativistic decoupling transformation. 

With gravitational coupling in a central, static field, 
the dominant attractive term is found to be described by 
the replacement £ —>• {§,1 — in Eqs. (HD) and (IMll . 
where £ is the energy operator defined in Eq. ([ 6 |) . This re¬ 
placement holds both for tardyons and tachyons and is a 
consequence of the structure of the Dirac-Schwarzschild 
Hamiltonian given in Eqs. (EH) and (EH- Namely, the 
dominant interaction in the high-energy limit is given 
by the anticommutator correction 5 {d?-p, (l — ^)} in 
the original Hamiltonians (before ultrarelativistic decou¬ 
pling) given in Eqs. (l22ll and (l42l) . The somewhat sur¬ 
prising observation that high-energy tachyons are at¬ 
tracted by gravitational fields finds a natural explanation 
in the energy-mass equivalence, and in the observation 
that both tachyons as well as tardyons travel at speeds 
very close to the speed of light in the high-energy limit. 
Indeed, the only difference in the effective high-energy 
Hamiltonians m and (1541) lies in the sign of the ki¬ 
netic term ±m^/( 2 £)), which is due to the changes in the 
dispersion relation for tardyons as opposed to tachyons. 
Higher-order corrections to the gravitational coupling are 
discussed in Appendix iDl 

The ultrarelativistic decoupling transformation should 
find applications beyond the description of gravita¬ 
tional interactions, for highly relativistic particles sub¬ 
ject to electromagnetic fields, and further applications to 
“nearly massless” electrons in graphene can be imagined 
(here, the “speed of light” is replaced by the Fermi veloc¬ 
ity VF, and dislocation potentials are added “by hand”, 
see Ref. [3^1. 


Appendix A: Dirac and Weyl Basis 

In the Dirac basis, we have 

12x2 0 


0 a 
a 0 


/3 = 


0 —l2x2 

0 


The 7 matrices in the Dirac basis are 7 ^^ = /3 and 


7 = 


/ 0 

5 1 

( 0 

l2x2 \ 

( 0 ) ’ 

7 = 1 

V 12X2 

0 j 


(Al) 


(A2) 







7 


In the Weyl basis, we have 

“w = ( o'" ^ ) ’ ^ o"" ) • (A3) 

We define the vector of dw matrices so that the “upper” 
solution includes the left-handed neutrino, whereas the 
“lower” spinor contains the right-handed Dirac antineu¬ 
trino [see Eq. (fT^ ]. The 7 matrices in the Weyl basis are 
7w = /3w and 



/ 0 

5 1 

f '—12x2 0 \ 

7w = 

( -a 0 ) ’ 

7w = 

( 0 12X2 j 


so that 7 ^ = —p. We notice that dw = Pw 7 w and also 
dw = —/3 E where 


E = 




(AS) 


The vector of Pauli spin matrices is denoted as a. Note 
that some authors define the vector 7 w with the opposite 
sign, which also reverses the sign of 7 ^ = i 7 ^ 7 w 7 w 7 w- 
Incidentally, the Coulomb coupling identifies particles 
(which are attracted) and antiparticles (which are re¬ 
pulsed). It is interesting to verify whether the interpre¬ 
tation is preserved under the transformation to the Weyl 
representation. We start from the Hamiltonian 

He = a- p- — , (A 6 ) 

r 


which describes a massless particles in a Coulomb field 
(here, Z is the nuclear charge number, while a is the fine- 
structure constant). Transformation to the Weyl repre¬ 
sentation is accomplished by the rotation 


He = UHcU-^ = -PE-p- 


Za 


—a ■ p — 

0 


Za 


0 


a ■ p - 


Za 


(A7) 


A comparison with Eq. (HH) reveals that states with pos¬ 
itive unperturbed energy (positive eigenvalue of the op¬ 
erator —<7-p for the upper spinor and positive eigenvalue 
of tr-p for the lower spinor) are attracted by the Coulomb 
field. By contrast, states with negative unperturbed en¬ 
ergy (negative eigenvalue of the operator -a-p for the 
upper spinor and negative eigenvalue oi a-p for the lower 
spinor) are repulsed by the Coulomb field. 


Appendix B: Operators 

We wish to explore the application of the operator 
I/f = —1/(E • p) to a reference state wave function. To 
this end, we assume that / = /(r) is a test function, and 


we defined the Fourier transform H and Fourier back- 
transform as follows. 


{TJ){k)= y d^rexp(-ifc-r)/(f), (BI) 

/ r|3p 

exp(ifc • r) f{k). (B2) 

One first multiplies the operator 

1/f ^-I/(S.fc) = -^ (B3) 

in Fourier space and then transforms back to coordinate 
space. 



r 


-1 


E- fc 
k"^ 


(-^/) (k) 


(r) 


d^k 

(^ 


d3./ 


r ^ ' k ^ik {f 

P 


/(r ). 
(B4) 


For a reference state with a special value kg of the wave 
vector [see Eq. mi 

fs{P) = MP\k^k.,t=o = ( ’ 

(B5) 

one has 

{J^fs){k) = ) S^^\k-P) (B 6 ) 

and so 

- - 1 .) . (B7) 

whose Fourier backtransform is 

(5^') “ ii ( “'P ) (*' I ■ 

This corresponds to the naive result that we obtain when 
interpreting the £ operator as an energy operator and 
applying it to the eigenstates of the free Hamiltonian, 
given in Eq. (na). 


Appendix C: Formalism for Gravitational Coupling 

We here follow the conventions used in Refs. i i 
for the flat-space and curved-space Dirac gamma matri¬ 
ces. Specifically, the flat-space and curved-space Dirac 
gamma matrices are distinguished in this Appendix using 
a tilde ( 7 ) and an overline ( 7 ) respectively. We draw in¬ 
spiration from the book and denote indices related to 
a local Lorentz frame (“anholonomic basis”) with capital 










Latin indices A,B,C,... = 0,1, 2,3. The curved-space 
Dirac gamma matrices 7 ^ satisfy the condition that 

{r(x),r(^)} = 2 r"(a;), (ci) 

where g^^{x) is the curved-space-time metric. The ^^{x) 
are expressed in terms of the flat-space Dirac 7 matrices 
7 "^ as follows, 

7^(x)=e^7^, lt,{x)=e^jA, (C2) 

where the are the coefficients which relate the locally 
flat Lorentz frame to the global space-time coordinates 
(the “vierbein”). Greek indices = 0,1, 2, 3 de¬ 

note the global coordinates. Latin indices starting with 
= 1,2,3,... are reserved for “spatial” global 
coordinates, which leaves /, J, iL, • • • = 1, 2, 3,... for spa¬ 
tial coordinates in the anholonomic basis. This nota¬ 
tion addresses some ambiguities which could otherwise 
result from other approaches [13, For exam¬ 

ple, unless the Dirac matrices are distinguished by over¬ 
lining or tildes, the expression 7 ^ could be associated 
with a flat-space matrix 7 ^^^ or with a curved-space 
matrix 7 *^^. We use the “West-Coast” convention for 
the flat-space metric, which we denote as = ijab = 
diag(l, —1, —1, —1). The curved-space metric is recov¬ 
ered as 

VAB = ^{ia,1b} , (C3) 

= ^{7f,{x),j^{x)} = e^e^r]AB, (C4) 

r-'ix) = \{Y{x),r{x)} = . (C5) 

For the curved-space metric around a gravitational cen¬ 
ter, we use the isotropic Schwarzschild metric in the Ed¬ 
dington reparameterization (43 |. i.e. 

9f,u = dmg{w^,-v^, -u^), (C 6 ) 

r" = diag(u;-2, -u-2, -u-2), (C7) 

u; = -^, ^=(1 + ^)'. (C 8 ) 

4r 


for a tachyon in curved spacetime as 

5 = y d'^x \/— det g{x) (dl) 

X i’{x)l^{x) Q7^(a;)Vp-7®(x)m^ , 

where tl’(x) 7^(x) takes the role of the “tachyonic adjoint” 
(see Ref. [Ig) and 

Vp = dp- rp(x), (C12) 

rp(x) = J uj^^ix) ^ab , dAB = ^ [ia,1b] , 

(CIS) 

{x) = V. = ef: d. FC^^ . (CM) 

Here, the F^^,^ are the Christoffel symbols, and the 
u!^^{x) are known as the Ricci rotation coefficients. 
Under a spinor Lorentz transformation with generators 
U^^(x), 

tjj'ix') = S{A{x)) ip{x) = exp aAB^ tp{x) 

(C15) 

we have covariance, i.e., Vltlj'{x) = VJ, [5'(A(x))'!/'(a;)] = 
{dl, - F(^) [S'(A(x))'0(x)] = 5(A(x)) Vj.'0(x), where 

the tranformed Ricci rotation coefficients F(, = 

j w'p (x) oab are calculated with respect to the trans¬ 
formed local coordinates. 

The curved-space fifth current 7 ^(a:) needs to be clari¬ 
fied. Adopting Eq. (18) of Ref. [13] for West-Coast sign 
conventions, one finds 

= 4^==^7'"(a:)7''(a()7^(a()7^(a(), (C16) 
4! y'-detg(a;) 

where e is the fully antisymmetric Levi-Civita tensor, 
with eoi23 = 1. We recall that the flat-space 7^ is 

=lb = j^ ^ABCDl^l^l^l^ ■ (C17) 

For the Schwarzschild geometry, one easily hnds 

Aet'g{x) = — , \/ — det 'g{x) = wv^ , (C18) 

7’5(x) = 75(0;) = 170717273 = 1^ = 1^ , (C19) 


For the Schwarzschild geometry, the vierbein coefficients j.e., the flat- and curved-space 7 ^ matrices are identical, 
read as follows. Variation of the action (ICllI) gives us 


II 

0 

Sp =SpV, 

(C9) 

— ifL 
w 

II 

(CIO) 


Here, 5^ = denotes the Kronecker-(5 (which is of 
course equal to unity for the two indices being equal and 
zero otherwise). 

In full analogy with the case of a “normal” massive 
Dirac particle (see Refs. [3,0); we write the Dirac action 


(iy'^ Vp - 7 ® m) i/f = 0 , (C20) 

which can be rewritten as 


i( 7 °)^ doip = ( 7 ° 7 V* + 17 ° y'" Fp -I- 7 ^m) . (C21) 

An explicit calculation of the Ricci rotation coefficients 
show that 


a ■ Vw 


d ■ Vu 

9 ’ 

WV^ 


(C22) 


7°7dp = 


2'w^v 
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which is in agreement with Ref. Si , where the a = 
7 ° 7 matrices are flat-space matrices. With the help of 
Eq. (IC21I) . one then finds that idtip = Hip, where 


w ^ ^ a - [p,': 

H = -a-p-\ -— 

V Iv 


[p^v]) 


■ w , 


(C23) 


needs to take into account the fact that in higher orders 
of the gravitational coupling constant, we cannot use the 
starting Hamiltonians as defined in Eqs. (EH) and (HI]). 
Instead, we must must use higher order approximations 
to the gravitational terms, which are otherwise neglected 
in Eq. (IC25I) (see also Refs. Si)- These lead to the 
initial Hamiltonians 


and /3 = 7 °. We now stretch the spatial coordinates with 
the he^ of the operator in analogy to the tardyonic 
case [J, Q , and find the 7 ®-Hermitian Hamiltonian 

H' = H + ™ w , (C24) 

2 t V i 

with H' = 7 ® 7 ^. Approximating w and v, according 

to Eq. (IC 8 D . to the first order in gravity, 

zc^l-g, ^-1 + g;, (C25) 

one finds the tachyonic gravitationally (TG) coupled 
Hamiltonian 

i?TG = ^{a-:P, (l-^)}+/375m (l - ^) • (C26) 

In the main body of the article, standard notation is ex¬ 
clusively used for the flat-space Dirac matrices (no over¬ 
lining and no tildes), i.e., we denote the 7 ^ as 7 ^^. 


Appendix D: Higher-Order Terms 

As discussed in Secs. El and ED the only differ¬ 
ence between the effective high-energy Hamiltonians for 
tardyons [Eq. (I41|l ] and tachyons [Eq. (15411 ] , derived in the 
main body of this work, is due the different dispersion re¬ 
lations for (free) tardyons and tachyons, while the gravi¬ 
tational interaction terms are identical to first order in 
and first order in l/£. This observation can be traced to 
the fact that the terms multiplying the kinetic operator 
and the mass in Eqs. (IMl) and (1421) . namely, X = I — 
and T = 1 — fulfill the relationship Y'^jX = 

One then easily reveals the cancellation mechanism for 
the terms of first order in Vs by treating X and Y in 
the non-transformed Hamiltonians (l22|l and (l42|l as con¬ 
stants. However, this does not imply that gravitational 
effects are the same for tardyons and tachyons, in higher 
orders of G (higher orders of r^). 

The ultrarelativistic decoupling transformation, keep¬ 
ing terms second order in rg, and to the first order in I jS 
(see AppendixEI, is expected to lead to differences in the 
gravitational interaction terms. In the calculation, one 


= 2 “ 




■ j3 m 


1 _ !± + Za. 

2 r 8 r 2 


for tardyons and 




^ , rs 9 r- 

■p,l - 

r 


16^2 j 




(DI) 


(D2) 


for tachyons. We then transform these Hamiltonians into 
the Weyl basis using the transform U defined in Eq. ®. 
Calculations become tedious and lengthy. One observa¬ 
tion in generalizing the decoupling transformation is that 
given a function / = /(r), then to first order of l/£ one 
finds 


^f£ + £fi=2f + ^[£,[£,f]]^^2f. (D3) 

As discussed in Sec. IHI A[ this is due to the fact that the 
two operators £ act only on the function /(r), and not 
on the wavefunction, thus they do not give “dominating” 
energy terms. After three iterations of the transform (per 
Hamiltonian), one finds for tardyons 



while for tachyons 



The final two terms in these Hamiltonians have opposite 
signs, indicating a difference in the gravitational interac¬ 
tion for tachyons and tardyons. 
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